Dynamics and relaxation in spin nematics 
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We study dynamics and relaxation of elementary excitations (magnons) in the spin nematic 
(quadrupole ordered) phase of S — 1 magnets. We develop a general phenomenological theory 
of spin dynamics and relaxation for spin-1 systems. Results of the phenomenological approach are 
compared to those obtained by microscopic calculations for the specific S = 1 model with isotropic 
bilinear and biquadratic exchange interactions. This model exhibits a rich behavior depending on 
the ratio of bilinear and biquadratic exchange constants, including several points with an enhanced 
symmetry. It is shown that symmetry plays an important role in relaxation. Particularly, at the 
SU(3) ferromagnetic point the magnon damping F depends on its wavevector k as T cx k 4 , while a 
deviation from the high-symmetry point changes the behavior of the leading term to T cx k 2 . We 
point out a similarity between the behavior of magnon relaxation in spin nematics to that in an 
isotropic ferromagnet. 

PACS numbers: 76.20.+q, 75.10.Jm, 75.40.Gb, 75.30.Ds 



I. INTRODUCTION 

Exchange interaction between atomic spins in solids 
at low temperatures usually leads to magnetic order- 
ing, that entails spontaneous breaking of the time re- 
versal symmetry. A large body of results in physics of 
magnetism have been obtained in the framework of a 
phenomenological theory, which describes the state of a 
magnet at low temperature by a constant length magne- 
tization vector M proportional to the average value of 
the spin operator (S) (or, in the more general case of a 
magnet with n sublattices, by several sublattice magne- 



tization vectors M a , a — 1, 2, ... 
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This approach 



naturally leads to the macroscopic description of the spin 
dynamics by means of the Landau-Lifshitz equations for 
the magnetization vector (or sublattice magnetizations). 

The assumption that the magnitude of the magnetiza- 
tion vector is conserved, \M\ — const, is justified in the 
case of a well pronounced long-range magnetic ordering 
with \M\ being not much different from its saturation 
value. In that case longitudinal oscillations that change 
\M\ have much higher energies than transversal modes 
described by the Landau-Lifshitz equations, and thus the 
longitudinal modes can be safely neglected in the low- 
energy, long-wavelength limit. However, the longitudinal 
modes may come down and become comparable in en- 
ergy with the transversal ones close to phase boundaries, 
where thermal or quantum fluctuations lead to a strong 
"spin contraction" so that | (S") | S*. In that case there 
might be other low-energy degrees of freedom as well, 
for instance, those related to the so-called spin nematic 
ordering 4 -^ that is described by the deGennes tensor or- 
der parameter 



Q 



ab 



S(S + l] 



Sab 



1 



(S a S b + S b S a ) 



built from quadrupole averages. Although the latter ex- 
pression makes sense only for S > \, nematic order is 
possible in spin-i systems as wellj£~— in that case the 
S > | spin operators S a above are understood as com- 
posed from S = o spins belonging to different lattice 
sites. In a usual magnetically ordered state, spin nematic 
(quadrupolar) order Q ab is trivially present as a "slave" 
of the primary magnetic (dipolar) order M. In contrast 
to that, in a purely nematic state the average (sublattice) 
magnetization M is absent even at zero temperature, and 
the only order is characterized by nontrivial quadrupole 
averages. The time reversal symmetry remains unbroken 
in the spin nematic phase. 

The static and dynamic properties of spin nematics 
have attracted much interest of researchers during last 
two decades JLHtSL This interest has received a new boost 
recently, mainly in the context of spinor Bose-Einstein 
condensates (BECs). 28 Optically trapped ultracold gases 
provide a unique highly controllable environment opening 
an exciting route for simulating a wide range of strongly 
correlated systems including quantum magnets. Ongo- 
ing experiments 2 ^— have already managed to reach the 
regime suitable to study magnetic properties of such sys- 
tems. Spinor gases are especially interesting because 
their internal degrees of freedom result in a rich physics, 
providing an opportunity to study quantum magnets 
with strong non-Heisenberg interactions. 

Spin-1 gas represents the simplest bosonic spinor sys- 
tem. Depending on interparticle interaction a 36 : 37 deter- 
mined by the s-wave scattering lengths ao,2 for the colli- 
sion channels with the total spin and 2, spin-1 BEC has 
ferromagnetic ground state for ao > a,2 (as in 87 Rb, see 
Ref.l38[). and nematic (polar) ground state for 02 > ao (as 
in 23 Na, see Ref . l39i) . For spin-1 atoms loaded in an opti- 
cal lattice a variety of phases has been predicted] 20 ' 40 " — 
At odd fillings, in the Mott insulator regime (intersite 
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hopping much smaller than the on-site interaction) spin- 
1 bosons on a lattice are described by the effective spin 
model that includes, in addition to the usual Heisenberg 
exchange terms (Si • Sj), strong biquadratic exchange of 
the type (Si ■ Sj) 2 . In two- and three-dimensional sys- 
tems, this biquadratic exchange leads to a spin-ncmatic 
state for a 2 > eto, while the case ao = a 2 exhibits an en- 
larged SU(3) symmetry with a highly degenerate ground 
stated 

In this paper we develop a general phenomenological 
theory of spin dynamics and relaxation in isotropic spin- 
1 magnets, focusing on the properties of the nematic 
phase. We compare the results of the phenomenological 
approach to those obtained by microscopic calculations 
for the S = 1 lattice model with isotropic bilinear and 
biquadratic exchange interactions. Our goal is to demon- 
strate the role of enhanced symmetry in the relaxation. 
Particularly, we show that at the SU(3) ferromagnetic 
point the magnon damping T depends on its wavevector 
k as r cx fc 4 , while a deviation from the high-symmetry 
point changes this behavior to T cx k 2 . Our formalism re- 
veals parallels between the general equations describing 
spin dynamics and relaxation in spin nematics and those 
proposed earlier by one of the authors^ 9 - for the dynamics 
and relaxation of magnetization in ferromagnets. 

The structure of the paper is as follows: in Sect. |TT] 
we describe the phenomenological theory of spin dynam- 
ics and relaxation, based on the Onsager relations, as it 
has been done earlier for ferromagnets, 49 In Sect. IIV| we 
discuss the microscopic calculation of magnon relaxation 
for the S = 1 bilinear-biquadratic model, and compare 
the results of the phenomenological and microscopic ap- 
proaches. Finally, Sect. [V] contains a brief summary. 



follows: 

n a = z *aKib z bi (2) 

where A", a = 1, . . .8 are the Gell-Mann matrices that 
are hermitean and have the following properties: 

\ a X P = ~8af)t + (d aM + lf aPl )\\ 
2 

^ab^a'b' = ^ab'Sa'b ~ ^ab&a'b' ■ (3) 

Here the tensor of structure constants / Q( g 7 is totally an- 
tisymmetric with respect to the permutation of any pair 
of indices, while the other tensor d Q ^ 7 is totally symmet- 
ric under such operations, see the Appendix. 

The octet n transforms according to the adjoint rep- 
resentation of the SU(3) group under a general uni- 
tary (SU(3)) transformation of the triplet of basis spin-1 
states. For two octets n and n', one can define a scalar 
product (n ■ n'), vector crossproduct (n A n'), and sym- 
metric vector product (n * n') as follows: 52 

Ti n' = n a n' a , (n * n') a — d a ^npn'^, 

(n A n') a = f a fanpn' r (4) 

For any octet n, quantities 

l2( n ) = n ■ n. Isin) = n ■ (n* n) (5) 

remain invariant under SU(3) basis tranformations z 1— ¥ 
Uz, and satisfy the constraint 

h(nf > 3I 3 (n) 2 . (6) 



II. PHENOMENOLOGICAL THEORY 
OF SPIN DYNAMICS AND RELAXATION 
IN S = 1 MAGNETS 

A. Description of spin-1 states and order 
parameters 

To describe the dynamics of spin-1 system, it is con- 
venient to use the formalism of SU(3) coherent states*^ 
The most general pure spin-1 state \ip)j at a single given 
site j is a linear superposition of three basis states \a)j 
with Sj\a)j = <y\a)j, o — 0,±1. It is convenient to write 
down this state in a "cartesian" basis of states 

IV") j = w 

a=l,2,3 

where | tl > = -J=(| - 1} - | + 1)), \t 2 ) = -^(| - 1} + 1 + 1)), 

1*3) = then the object Zj = (zji,Zja,Zj3) transforms 
as a vector under usual (SU(2)) rotations. 

Instead of the three-component complex vector z, 
one can use a different spin-1 state parametrizatio n 50 i 51 
through the eight-component real vector n, defined as 



The density matrix p of a single spin can be expressed 
through the octet n: 

(ta\p\t b ) = ls ab + ^n a X^. (7) 

For a normalized pure state, the condition p 2 — p trans- 
lates into the following constraints on n: 

4 2 
n 2 = —, (n * n) — —n. (8) 

i.e., I2(n) — 4/3 and h(n) = 8/9, so (JSJ) becomes an 
equality. One can show that the constraints © reduce 
the dimension of the n-space to four i 12 ' 50 

For a mixed state, Eq. ([7]) retains sense as it is the most 
general expression for a 3 x 3 matrix with unit trace. The 
constraints © generally do not hold for a mixed state. 
However, under a unitary rotation p(n) 1— > U'f>(n)U = 
p(n') the transformation n H> n' is still determined by 
the adjoint representation of SU(3), so l2,s(ii) = h,3(n') 
remain invariant, but are no more fixed at the pure state 
values ([S]). The dimension of the n-space in the general 
case of a mixed state (J7J is equal to sixi 12 ' 50 
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The components of the octet n correspond to the fol- 
lowing on-site spin averages: 

n 2 = {S z ), n B = -(&>), n 7 = (S x ), 

m = (S x S y + S y S x ), n 4 = -(S X S Z + s z s x ), 

n 6 = (S V S Z + S Z S V ), (9) 

n 3 = ((S x ) 2 - (S v ) 2 ), n 8 = V3(((S Z ) 2 ) - 2/3), 

which can be split into the magnetization (dipole) part 
m and nematic (quadrupolar) one d, 

n = m + d, 

m = (n 7 , -n 5 ,n 2 ), d = (m, n 3 , n 4 , n 6 , n 8 ). (10) 

It is easy to see that under time reversal operation T, 
the magnetization m changes sign, while the quadrupolar 
part d remains invariant: 



Tm 



-m. Td = d. 



(11) 



B. Equations of motion for pure spin S = 1 states 

The effective Lagrangian of a spin-1 system on a lat- 
tice, expressed in terms of the complex unit vector z, 
takes the form 

C = Y^i(z*j ■ d tZj ) - W({z*, Zj }), (12) 

3 

where is the energy that can be also expressed through 
the octets {rij}. We will assume that one can pass to the 
continuum description just by declaring z (or, alterna- 
tively, n) a smooth held. Then the Lagrange equations 
of motion for z* , z take the form 



(13) 



.dz* _ SW 
l ~dT ~ ~~Sz~' 

where the energy W[z* , z] = W[n] is now a functional of 
either the complex vector held z or the real octet held n, 
and S/Sz denotes a variational derivative. Introducing 
generalized fields 



H 



SW 
Sn ' 



h = - 



SW 
~Sz~' 



(14) 



one can write 



dn a .. . 
— = -i{h a \ ab z b -c.c), 

which finally leads to the equation of motion for n: 



dn 
~dt 



2(nAH). 



(15) 



It is worth noting that if W depends only on magnetiza- 
tion m, (|15[) reduces to the well-known Landau-Lifshitz 
equation without dissipation, dm/dt = — (m x SW/Sm), 
which describes spin dynamics in ferromagnets. Thus, 
(|15p can be viewed as an extension of the Landau-Lifshitz 
equation that includes the dynamics of quadrupolar de- 
grees of freedom. 



General form of the equations of motion and 
relaxation terms 



General analysis of dissipative forces in spin-1 nematics 
can be carried out closely following the approach devel- 
oped for ferromagnets4^ Equations of motion can gener- 
ally be cast in the form of Onsager equations 



dn a 
~dt 



= A a p(n)Hp, 



(16) 



where the effective fields H, defined in (|T4|), play the role 
of generalized forces that arise at a deviation from equi- 
librium, and A a p are the kinetic coefficients that gener- 
ally shall be understood as operators acting on H (phys- 
ically this corresponds to taking into account the spatial 
dispersion). 

Since different components of the octet vector n have 
different properties under time reversal, the Onsager 
theorem about the symmetry of the kinetic coefficients 
A a /3 (n) gets slightly modified. 53 If one introduces the 
factors e a = ±1 describing the symmetry of the a-th 
component of n under time reversal, 



(17) 



then the kinetic coefficients satisfy the following reci- 
procity relations: 



Ap a (n) = e a epA a p(Tri). 



(18) 

It is easy to see that only the symmetric part A L 

of the tensor A contributes to relaxation. Indeed, the 
dissipative function Q can be written as 



ldW _ 1 
^ ~ ~2~dT ~ 2 



dt 



1 



dx H a A a pHp 



The antisymmetric part A^ of the kinetic coefficients 
tensor makes no contribution into dissipation and thus 
corresponds to the purely dynamic part of the equations 
of motion. As we have already established, those equa- 
tions are given by (|15l) . and thus one obtains 



1 



dxHj^Hp. (19) 



(20) 



Under the assumption of weak dispersion, one can ex- 
pand the symmetric part of A in spatial derivatives, re- 
taining only up to quadratic terms. 

Assuming further that the crystal structure of the mag- 
net has an inversion center, one can ignore the term linear 
in derivatives, so finally we obtain: 



( n ) = ^W( n ) - KpM'{n) 



d 2 



dxi dxi> 



(21) 



where \ a p and X a p,w satisfy the reciprocity relations 
similar to (|18[) with respect to permutations of a and /?, 
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and \ a p,w are symmetric with respect to permutations 
of I and V . The dissipation function in this approxima- 
tion will contain terms quadratic in the effective fields 
and their derivatives, 



= dx 



Hp + X a 0,w 



dH a dH, 



dxi dxi 



5}' (22) 



and equations of motion including relaxation terms ac- 
quire the following general form: 



dn a 
~~dt 

Ret 



Zfap-fnpH^ + R a 



— X a pHp — \ a p,W 



6H C 



dxi dxy 



(23) 



III. PHENOMENOLOGY APPLIED TO 
BILINEAR-BIQUADRATIC MODEL 

We apply the general formalism presented above to the 
bilinear-biquadratic model described by the Hamiltonian 



u 



(ij) 



(ij) 



(24) 



where Si are spin-1 operators at the i-th lattice site, and 
the sums are over nearest neighbors. For simplicity, we 
assume that the lattice is cubic, and the lattice constant 
is set to unity. 

The above model describes the most general isotropic 
(SU(2) invariant) exchange interaction between two S — 
1 spins. At Ji = J2 the symmetry of the model is en- 
hanced to SU(3). We consider the nematic region in the 
vicinity of the ferromagnetic SU (3) point, so it is conve- 
nient to set 

Ji = J(l-5), J 2 = J, J>0, < <5 < 1/2. (25) 

Introducing the continuum field n = m + d, one can 
write down the energy W = (H) as 

W = J d 3 x{(Z/2)[-^Jn 2 + 5Jm 2 ] 

+ ^[(Vd) 2 + (l-25)(Vm) 2 ]}, (26) 

where Z = 6 is the lattice coordination number, and 
(Vn) 2 = *Yln{dn/ dxi) 2 . For pure states, the constraints 
© fix the length of n; thus, for 5 > 0, in the ground 
state the magnetization m vanishes, which corresponds 
to a uniaxial spin nematic. 

At finite temperature T, expression (l26t has to be re- 
placed by the free energy, T — W — TS, where S is the 
entropy of the system. Below we use a phcnomcnolog- 
ical expression for the free energy written in the spirit 
of Landau's theory of phase transitions in a form of ex- 
pansion on powers of order parameters. This expression 
generally can be constructed from all possible invariants 
of the corresponding symmetry group. 



We limit ourselves to the case S < 1/2, when local 
spin correlations are ferromagnetic. For 6 > 1/2 the 
nearest-neighbor spin-spin correlations change their char- 
acter to antiferromagnetic, and the effective continuum 
theory has to be modified (n has to be split into a uniform 
and staggered components that become smooth fields in 
the continuum theory^) . The detailed discussion of this 
interesting case is going far beyond the scope of this arti- 
cle, and below we will limit ourselves to the case S < 1/2, 
paying main attention to the vicinity of ferromagnetic 
SU(3) point, where 5 < 1. 

As a consequence of the assumed cubic symmetry of 
the lattice, the tensor of dissipative constants, that en- 
ters relaxation forces depending on derivatives, must be 
diagonal in its space indices, 

A a /3,2Z' = XapSw- 

Now let us focus separately on the highly symmetric case 
5 = and then study what happens at deviations from 
this point. 



A. 5(7(3) symmetric case 

At 6 = the system is S'J7(3)-symmetric, which dic- 
tates that at this point both tensors of dissipative con- 
stants are diagonal in the octet indices, \ a p = \8 a p and 
Xap — ^5 a p. The dissipative force R is thus determined 
by just two constants A and A: 



R=XH -XV 2 H. 



(27) 



The free energy J 7 will generally, in addition to the 
gradient term, contain some function / \l2,Iz) of the 
SU(3) invariants l2(n) and Ia(n), that will have a min- 
imum at certain equilibrium values of -^(^o) — ^2 



li 0) and 



^3(^0) = ^3°'- Keeping only the lowest order terms in 
the expansion of /(/2, ^3) around its minimum as well as 
in the gradient expansion, we obtain 



T= \ d 



(n 2 -I^) 2 , (n.(n*n)-/n 



(0h2 



4X2 



6X3 



(28) 



where \2 and X3 play the role of longitudinal suscepti- 
bilities. The effective field H = —5F/8n thus takes the 
form 



H = -[(n 2 -I^)/ X2 ] 



n 



(29) 



- [(n ■ (n *n) - jfVxa] (n * n) + ^V 2 n, 

where the first two terms do not contribute to the equa- 
tion of motion (| 15[) in absence of dissipation (note that 
nA(n*n)=0). 

Further, at 8 = the quantity N = J d 3 xn is 
conserved. As a consequence, the equation of motion 
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for n should have the form of a continuity equation 
dn a /dt + divII Q = and so the dissipative force R must 
have the form R a = dH a b/dxt- It is clear from (f27f and 
(|2l))) that R can be represented in such a form only if 
A = 0. Thus, at the SU(3) symmetric point S = we are 
left with a single relaxation constant X: 



R=-XV Z H. 



(30) 



The above argument is fully similar to the treatment of 
the exchange approximation in ferromagnets4^ 

Let us calculate the magnon damping in spin nematic 
at the SE/(3)-symmetric point. Without loss of general- 
ity, the uniform ground state can be chosen as n = no 
with 



n = no(0, 0,0, 0,0, 0,0,-1). 



(31) 



This state corresponds to a uniaxial spin nematic with 
zero magnetization and the nematic director along the 
z-axis. To visualize it, one may think of an ellipsoid with 
the same symmetry as for the set of quadrupole averages 
(SiSj + SjSi). For the above state, such an image is 
a squeezed ellipsoid of rotation, with the ellipsoid axis 
being directed along the z axis, and with the thickness 
oc (S 2 ) smaller than the diameter cx (S 2 ) = (<S?). At 
zero temperature no = 2/V3 and thus (S 2 ) = 0, i.e., 
the ellipsoid degenerates into a zero thickness disk. Due 
to the high symmetry of the system, any other choice of 
no that can be obtained from ([31) by a general SU(3) 
rotation would give the same physical results. We have 
chosen the ground state ([31]) because it has the simplest 
form in the octet representation. 

We consider small deviations from the ground state, 
n = no + n, linearize equations of motion (|23|) in n, and 
look for eigenmodes in the form n = jje % ^ k ' 31 ~ ut > . Then 
the damping F& is obtained as the imaginary part of the 
frequency uik = &k — i^k- Since we are interested only 
in the case of weak inhomogeneities, one can expand flk 
and Tfe in powers of the wave vector k and retain only 
leading terms. 

It turns out that the dynamics of the components rji, 
r]2, T]3, and 778 is purely diffusive: the real parts flk of 
the corresponding frequencies vanish, and the damping 
behavior is given by 



ri = r 2 = r 3 = ^xjk 4 , 



(32) 



r* 



CXk 2 + -XJk 4 , C 



The above modes are longitudinal, they correspond to 
hydrodynamic relaxation of inhomogeneities in the dis- 
tribution of n and are not connected to any propagating 
excitations. Let us briefly discuss their physical meaning. 
The variable r\i = m z describes the z component of the 
magnetization, and two variables 771 , 773 correspond to 
quadrupolar averages built from projections of the spin 
in the (xy) plane. The dynamics of these variables can 



be viewed as a change of the spin length coupled with the 
rotation of the system around the z axis. Such a geom- 
etry corresponds to the longitudinal mode found in the 
same model at S < (i.e., in the ferromagnetic state), 
where this mode acquires a finite frequency and thus be- 
comes propagating^ The last diffusive mode r] S deter- 
mines the evolution of the variable (m 2 ) and describes 
the relaxation of the de Gennes order parameter to its 
equilibrium value. In the S'C/(3)-symmetric case there 
is only inhomogeneous relaxation (all damping constants 
contain k). As we shall see below, a deviation from this 
high symmetry point leads to the appearance of a homo- 
geneous relaxation for quadrupolar degrees of freedom, 
but not for the total magnetic moment. 

The rest of the modes are dynamical and correspond 
to propagating generalized magnons that unite dipolar 
(spin) and quadrupolar excitations. The dynamics of 774 
and «5 is coupled, and the same holds for the pair r]Q, r]7. 
The corresponding cigcnfrcquencies are 



U>( 4 



= ±n 



(45) 



iT 



(45), 



w (67) = ±0(67) _ £T(67)>(33) 



(45) 



(67) 



Jnov3 



k 2 . 



1 



(45) 



(67) 



XJk 4 . 



One can see that at small wavevectors the magnon damp- 
ing behaves as F^ cx k , while the real part of the fre- 
quency Qk cx k 2 3> Ffc, so magnons remain well-defined 
excitations. 



B. Spin nematic away from the 5(7(3) point 

Consider a system in the nematic phase (1/2 > 5 > 0). 
In the ground state the magnetization m vanishes, so the 
equilibrium value n = no = do lies completely within the 
quadrupolar subspace. The symmetry group is now re- 
duced to SU{2) which is a subgroup of SU(3). Under 
rotations, m and d transform according to D 1 and D 2 
representations of the SU (2) group, respectively. Thus, 
m 2 and d 2 are invariants that can separately enter the 
free energy. The SU(3) invariants ^(n) and ^(n), of 
course, remain invariant under any subgroup transforma- 
tion; since I2 = m 2 + d 2 , we shall only take into account 
ia(n) as another independent invariant. We will assume 
that the free energy has the form 



T = 



2 \2 



SJm 2 



(d 2 -n 2 ) 



4xd 



1 ■ (n * n) - h(n Q )] 2 
6X3 

[(Vd) 2 + (l-2^)(Vm) 2 ]} ! 



(34) 



where we have neglected the terms of higher than 
quadratic order in m and (d — no). 

Similar to ([T0|) . we can divide the effective field H 
into two components H m and Hd lying in the dipolar 
and quadrupolar subspaces of the octet space, and trans- 
forming according to D l and D 2 representations of the 
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SU(2) group, respectively. Then, the symmetry dictates 
the following structure of the dissipative constants ten- 
sors: 

R = ■ H m — A m V H m + XdHd — AdV 2 Hd, (35) 

where the vanishing constant in front of H m is the con- 
sequence of the fact that the total magnetization M = 
J d 3 x m is an integral of motion. From our analysis of 
the SU (3)-symmetric point it also follows that 

A d -> 0, A m -> \ d at 5 -> 0. (36) 

Eqs. fl35l36[) are valid on both sides of the SU (3) point 
(i.e., at any sign of S), while Ad is strictly non- negative. 
If one assumes Ad to behave analytically as the function 
of S, then it follows that 

A d = 0{8 2 ) at 5 -> 0, (37) 

however, on general grounds one cannot exclude the pos- 
sibility of singular behavior of Ad (6). We will check the 
above assumption and its consequence (1371) in the next 
section by comparing it with the results of microscopic 
analysis. 

Let us now calculate the magnon damping away from 
the SU(3) point. Without loss of generality, we assume 
the nematic ground state with the director along the z 
axis, i.e., no of the form ([3~Tj) . We linearize equations of 
motion in the deviation r\ = n — no and find the complex 
eigenfrequencies lj — fl — iT, as we have done before 
in the SU(3) case. We find that there are again four 
longitudinal, purely diffusive (ft = 0), decoupled modes 
Wit r /2, f]3i an d % j with the following linewidths: 

ri = r 3 = ijfc 2 (Ad + Adfc 2 ), 

T 2 = X m Jk 2 [ZS + ~ (1 - 2S)k 2 ] , (38) 

/ 1 \ ~ / 2 n 2 \ 

T 8 = (c+-Jk 2 )(\ d + \ d k 2 ), C = nl{— + 

v 2 I \\d Xz> 

Comparing the above expressions with their SU(3) coun- 
terparts in Eqs. (|32|). one can see that the main effect 
of breaking the SU (3) symmetry is the appearance of a 
homogeneous relaxation in the rj 8 quadrupolar order pa- 
rameter proportional to ((S z ) 2 — 2/3): k = fluctuations 
of (S 2 ) decay with the characteristic relaxation time 

T8 55 T^iho) = iki (39) 

which diverges at 6 — > according to Eq. (|37|) . 

As before, there is no homogeneous relaxation of the 
magnetization, because the SU(2) symmetry remains in- 
tact and the total magnetization is conserved. For the 
other three diffusive modes 771,2,3, there is a small (of 
the order of S) splitting between the damping constants 
of the quadrupolar modes (771,3) and the magnetization 
mode (772), and the leading term in the damping at small 



wave vectors behaves as k 2 (in contrast to the k 4 behavior 
in the SU(3) case). 

The remaining two degenerate eigenmodes correspond 
to coupled (774,775) and (776,777) oscillations, and describe 
propagating magnons. The corresponding eigenfrequen- 
cies ft us) = f2(67) and dampings Tu$) — ^(67) are given 
by 

0(45) = "(67) = (2Z6 + (1 - 2S)k 2 ) , 

Jk 2 r 

r ( 4 5) = r (67) - — {Ad + 2Z5X m (40) 
+ [Ad + (l-25)A m ]fc 2 }. 

It is easy to see that away from the SU(3) point magnons 
acquire linear dispersion, in agreement with previous re- 
sults obtained by different methods^^ii^ and the leading 
term in their damping at small wavevectors behaves as 
k 2 . At S — > 0, taking into account (f31)|). we recover our 
results obtained above for the SU(3) point. 

IV. MICROSCOPIC ANALYSIS OF THE 
BILINEAR-BIQUADRATIC MODEL 

Our goal is to compare the results for magnon damping 
in the bilinear-biquadratic model (|2"4")l , obtained above 
within the phenomenological approach, with a calcula- 
tion from first principles. We would like to pass to 
the second-quantization formalism for the description of 
magnons. It is convenient to break up the complex vector 
z, describing the coherent state ([I}, into two real vectors 
representing its real and imaginary parts z = u + iv. 
When the normalization condition and the arbitrariness 
of the phase factor are taken into account, these vectors 
satisfy the conditions: 

u 2 +v 2 = l, u-v = 0. (41) 

In terms of the variable vectors u and v the Lagrangian 
of the system can be written as 

C = ^hJ^VjdtUj - W, 

j 

W = ^ {2(^2 - Jl) [(UiUj)(ViVj) - (UiVj)(ViUj)} 

-y [i^iUj + ViVjf + (iiiVj - ViUjf] }. (42) 

The spin nematic state, favored at J 2 > J\ > 0, corre- 
sponds to the following condition: vectors Uj at all the 
sites are parallel, Uj = Uo, \uq\ = 1, and Vj =0. (There 
is also an alternative configuration obtained by the sub- 
stitutions uj i-> Vj,Vj h> — Uj, but those two states are 
physically identical). 

Consider small deviations of the variables u and v from 
the spin nematic state. Assume, for the sake of definitc- 
ness, that Uq = e z . Then, according to (T4"Tj) . v z and 
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u z are dependent variables quadratic in the remaining 
components u x>y and v x y . We can thus expand the La- 
grangian in powers of the independent variables u x ,y and 
v x ,yi keeping terms up to the quartic order. From the 
kinetic part of the Lagrangian it is obvious that u x and 
u y can be chosen as the coordinates, while —2hv x and 
—2hv y play the role of the corresponding canonical mo- 
menta. It is easy then to pass from the Lagrangian to 
the Hamiltonian. Passing from variables Uj, Vj to their 
Fourier amplitudes u kl v k and performing quantization, 
we can write the standard representation of the coordi- 
nates and momenta in terms of the Bose creation and 
annihilation operators: 



^x,k 



2A t 



-(a\ + o_h), U X) fe = i\ 



Uy ' k = v^r^** + b ~ k ^ Vy ' k = *vW^ tfc ~ b - k) ^ 

where A k = 2h^/a k / (3k, and we use the notation 



a k = ZJ 2 [1 - 7 (fe)] , k = J 2 kfa(k) 



k 2 = 2Z(l - Jx/Ja) = 7 (fc) 



z ^ 



Oik, 



(44) 



£ being the set of vectors connecting a site of the lattice 
to its nearest neighbors. 

The Hamiltonian decomposes into a sum of quadratic 
and quartic terms: H = H 2 + H4. The quadratic part 
H 2 becomes diagonal, 



#2 = ^ £k{a k a k + b k b k ), 



and the magnon dispersion relation reads 



£k = V otkfik- 



(45) 



(46) 



The spectrum contains two degenerate magnon modes 
with orthogonal polarizations, that are gapless at k — > 

I 



in accordance to the Goldstone theorem. For J\ > J2 
this spectrum becomes unstable, which corresponds to a 
transition to the ferromagnetic state. 

In this paper, we are interested in the low-energy 
dynamics. In the vicinity of the ferronematic SU (3) 
point Ji = J 2 (i.e., (5 = 0) this is translated into the 
long- wavelength approximation k <C 1, so we obtain 
7(fc) wl- k 2 /Z and 



£fc 



J 2 k\/k 2 , 



k 2 at k < 1. 



(47) 



The quartic Hamiltonian H4 is rather cumbersome, 
but for our purpose of calculating the magnon damping 
in the lowest order of the perturbation theory it can be 
simplified substantially. Indeed, the analysis shows that 
decay processes (one magnon decaying into three) gives 
a small contribution at high temperatures, see the note 
after Eq.(|52|l below. Thus, the main contribution to the 
damping comes from terms containing an equal number 
of creation and annihilation operators. When only such 
terms are taken into account, the Hamiltonian takes the 
form 



1 

N 



E 

1,2,3,4 



Al+2-3-4 

+ F(a\a\a 3 a4 : + b\b\bsbi) + h.c. j 



(48) 



where for the sake of brevity we use the shorthand nota- 
tion 1 = fci, etc. 

The amplitudes "J, and F in Eq. (|4"5)l depend on 
all four magnon momenta. 4> corresponds to processes of 
conversion between pairs of magnons with different polar- 
izations, while the other two amplitudes describe magnon 
scattering. Further, in our case (in the vicinity of the fer- 
ronematic SU(3) point) for low-energy processes the mo- 
menta of all the magnons participating in the process are 
small, and one can use the long-wavelength asymptotic 
expressions for the corresponding amplitudes: 



$1 



4 



234 



:{ - k 2 (4kik 2 k 3 k 4 + nt, 2 ^t,i) + ( fclfc 3 + ^^(Vi^a - ^1,2^3,4) 



1234 



32- v /£i£ 2 £ 3 £4 

+(feife 4 + k 2 k 3 )(ri^ 2 r]- A + ^2^,4)} 

■■{k 2 ((ii 3 fi2A - 4fcifc 2 fc 3 fc4) - {kik 2 + fe 3 fc4)(?7^3?7^ 



•'2 



16 v /£i£ 2 £3£4 

+ (kiki + k 2 k 3 )(r]+ 3 r]+ i - nf,3f4,S\ 



^1,3^2,4) 



-F1234 — 



J 2 



Ml,3M2,4 

\ kli^kxk^kaki + kik^vti + kshvt 

32 v /£i£ 2 £3£4 

2(fcifc 2 + fo^Ml^M^ - 2 (klk 3 + k 2 ki) fJ-i^^A 

I 



^1,2/^3,4) + 2 ^t.2%,4 + Ml, 2/^3,4) 



(49) 



Here we have used the shorthand notation 



The magnon damping T is determined as the imaginary 



'^nvm-j Mn.rn — k n 8m i ^n&r) 

$n = yj k,Q ~\~ h = -\- k^ ~h ^3 ~h k^. 



(50) 
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part of the self-energy calculated in the second order of 
perturbation theory. The expression for the damping can 
be written in the following form: 

r^ = ^sinh(§)^($ 2 + 2^ + 4F 2 ) 
p i 

S (gfe + £p - £q - gfc+p-q) 

sinh (||) sinh (|f,) sinh (^ff^) ' l ° ' 

where T is the temperature in energy units, and the ar- 
guments of the amplitudes (|^|) are taken as k\ = fc, 
k-2 = p, &3 = q, &4 = k + p — q). We will consider the 
case of high temperatures, 

e k < T < J 2 , (52) 

then the main contribution to the damping comes from 
thermally excited magnons with wave vectors p much 
higher than the magnon wave vector k. This justifies ne- 
glecting magnon decays in the interaction Hamiltonian: 
although such processes are allowed by the energy and 
momentum conservation, their phase volume is propor- 
tional to k 2 and thus their contribution is small. 

We will also assume that we are not too far away from 
the SU{2>) point, so that 

T> J 2 kl = 2Z{J 2 - J x ). (53) 

At the SU{2>) point (fco = 0, 5 n — k n ) the amplitude 
$ vanishes, and in the immediate vicinity of the SU(3) 
point determined by (|53[) it is much smaller than the 
other two amplitudes, so <fr can be neglected in the cal- 
culation of the damping. 

Then, in the leading order in (J 2 — J\) and k/p, one 
can drop ko in the dispersions of thermal magnons, and 
the amplitudes in (j5"Tj) can be replaced by 

j3/2 

1-1234 = 2^1234 = -J^ih+SJikx-kz+ks-kl), (54) 

where we have dropped contributions proportional to 
k\jk\k 2 <C 1. Further, if k ^> fco- then fco can be ne- 
glected in the equation for the mass surface as well, and 
then the result for the damping under the assumptions 
(f52")l . ([53"]) can be easily obtained as follows: 

I\t = J -jC KT k 2 (2k 2 + k 2 + 0(A#)) at fc 2 «fc 2 , 

(55) 

where Ck,r is a dimensionless factor that weakly (loga- 
rithmically) depends on the wave vector k: 

^=8^{( ln S4) 2 J- (1 - ln2 ) 2 }- < B6 > 

At S = this result reduces to the one found before 
for the SU (3) symmetric caseM The appearance of the 
logarithm above is similar to the well known result for 
isotropic ferromagnets. 56 



Comparing (|55"|) to the phenomenological result (|40l) . 
one can see that it is possible to match those two expres- 
sions if one sets the correspondence 

X m ps X d h+ C k ,T, \d = ((J 2 - Ji) 2 ) , 

l m -Xd = 0(J 2 - Ji). (57) 

We see that our conjecture (l37l) . which stems from the as- 
sumption that Xd is an analytical function of ( J 2 — J\), is 
confirmed by the results of the microscopic calculations. 

V. SUMMARY 

To summarize, we have studied the properties of ele- 
mentary excitations in the so-called spin nematic phase 
of spin-1 systems that is characterized by nontrivial 
quadrupole order in the absence of local spin averages. 
We have developed a general phenomenological theory of 
spin dynamics and relaxation for spin-1 systems, which is 
based on the equations of motion for the eight-component 
real vector n uniting the magnetic (dipolar) and nematic 
(quadrupolar) order parameters. Our approach is well 
suited to emphasize the role of enhanced symmetry in 
the relaxation and is in spirit similar to the phenomenol- 
ogy describing relaxation in ferromagnets 4^ 

The developed theory has been applied to the spe- 
cific S = 1 lattice model with isotropic bilinear and bi- 
quadratic exchange interactions, which is relevant for the 
physics of ultracold spin-1 Bose gases in optical lattices. 
In the space of the model parameters, there is a spe- 
cial point that exhibits an enhanced symmetry, namely, 
the usual SU(2) (rotational) symmetry is enhanced to 
SU(3), and the immediate vicinity of the SU(3) point 
is occupied by the spin nematic phase. We show that 
the behavior of the leading term in the dependence of 
the magnon damping on wavevector k changes from k A 
to k 2 as one moves away from the SU(3) point and the 
symmetry is lowered to SU(2). Those predictions of the 
phenomenological theory are consistent with the results 
of microscopic calculations. We also show that break- 
ing the SU(3) symmetry leads to the appearance of a 
homogeneous (k — ¥ 0) relaxation in the diffusive (non- 
propagating) quadrupolar mode that describes fluctua- 
tions of (S 2 ). 
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Appendix: SU(3) algebra and properties of octet 
vector products 

For the sake of the reader's convenience, we list here 
the basic algebraic relation o 52 i 57 between the Gell-Mann 
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matrices X a and the structure tensors f a fs~i, d a p-y 



d a py defined by their nonzero components as follows: 



The generators of the SU(3) group can be chosen in 
the explicit matrix representation known as Gell-Mann's 
matrices: 



/l23 — 1) /458 — fe78 — , 
/l47 = /246 = /257 = /345 = /5I6 = /637 = 77 



(A.2) 



1 

2' 



Ai 



A 4 






(0 


-i 0\ 






i 





|.*-( 




v° 


0) 






(0 


-i\ 
















\> 


I 







(A.l) 



The algebraic properties of Gell-Mann's matrices are 
given by Eq. ([3]) , with the totally antisymmetric structure 
tensor / Q( 3 7 and the totally symmetric structure tensor 





— <^228 


— ^338 


— — d 8 88 : 


G?448 


= ^558 


= ^668 


= ^778 = 


^146 


= ^157 


= <^256 


= ^344 = 




= d 2 47 


= ^366 


= ^377 = 



1 

1 



2>/3' 



(A.3) 



^355 
1 



With the help of the above structure constants, one can 
define the symmetric and antisymmetric vector products 
as in Eq. @. 

If one splits the full space of octet vectors into the 
magnetization subspace V m and the quadrupolar sub- 
space Vd, see Eq. (fT0|) . then it is straightforward to check 
that the vector products have the properties 

(dAd)eV m , (dAm)eV d , (m A m) £ V m , 
(d*d) e V d , (d*m)£ V m , (m*m) £ V d , (A.4) 

where m £ V m and d £Vd- 
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